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ABSTRACT 

We analyze the correlations between central dark matter (DM) content of early-type galaxies and 
their sizes and ages, using a sample of intermediate-redshift (z ~ 0.2) gravitational lenses from the 
SLACS survey, and by comparing them to a larger sample of z ~ galaxies. We decompose the 
deprojected galaxy masses into DM and stellar components using combinations of strong lensing, 
stellar dynamics, and stellar populations modeling. For a given stellar mass, we find that for galaxies 
with larger sizes, the DM fraction increases and the mean DM density decreases, consistently with 
the cuspy halos expected in cosmological formation scenarios. The DM fraction also decreases with 
stellar age, which can be partially explained by the inverse correlation between size and age. The 
residual trend may point to systematic dependencies on formation epoch of halo contraction or stellar 
initial mass functions. These results are in agreement with recent findings based on local galaxies by 
Napolitano, Romanowsky & Tortora (2010) and suggest negligible evidence of galaxy evolution over 
the last ~ 2.5 Gyr other than passive stellar aging. 

Subject headings: dark matter - gravitational lensing - galaxies : evolution - galaxies : galaxies : 
general - galaxies : elliptical and lenticular, cD. 



1. INTRODUCTION 

Early-type galaxies (ETGs) are among the most mas- 
sive systems in the Universe. They are on average metal- 
rich, dust and gas poor, and formed their stars in early, 
rapid events (e.g., Thomas et al. 2005; Fontanot et al. 
2009). The lack of easily interpretable dynamical trac- 
ers, such as the cold gas in spiral galaxies, has made 
the mass mapping of these systems very difficult in the 
outer regions where dark matter (DM) is expected to 
be dominant, although discrete tracers (planetary neb- 
ulae and globular clusters; e.g., Napolitano et al. 2009; 
Romanowsky et al. 2009) are providing a clearer view 
of the DM halos. The mass content of these galaxies' 
central regions has on the other hand been extensively 
investigated (e.g., Gerhard et al. 2001; Cappellari et al. 
2006; Tortora et al. 2009, T+09 hereafter), with build- 
ing evidence that the central DM fraction (/dm) is an 
increasing function of total stellar mass (A/*), providing 
the main driver for the tilt of the fundamental plane (e.g., 
T+09). 

New insights into the galaxy assembly process are 
emerging from joint analyses of the structural and star 
formation properties of nearby ETGs (e.g., Gargiulo et 
al. 2009; Graves & Faber 2010; Napolitano, Romanowsky 
& Tortora 2010, NRT10 hereafter). A key discovery of 
NRT10 is an anti-correlation between /dm and galaxy 
stellar age, such that older galaxies (at a fixed M*) have 
lower /dm- In the context of ACDM halos, this trend 
can be partially explained by an anti-correlation between 
galaxy sizes and ages, with the remaining effect appar- 
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ently driven by variations in star formation efficiency, 
stellar initial mass function (IMF), or DM distribution 
(e.g., adiabatic contraction, AC hereafter). As discussed 
in NRT10, such correlations would have deep implica- 
tions for the assembly histories of ETGs, but critically 
need to be confirmed by independent analyses. 

Gravitational lenses offer a unique tool to map the 
mass profile in galaxies over a range of redshifts. The 
database of lenses is growing quickly thanks to ongoing 
surveys (e.g., SLACS: Auger et al. 2009, A+09 hereafter; 
COSMOS: Faure et al. 2008). In particular, the SLACS 
survey has collected >80 secure lenses, which is a sample 
comparable to the total number of lenses discovered since 
the late 70s from other campaigns or by serendipity (e.g. 
Covone et al. 2009). Here, we will use the SLACS sample 
data to extend the analysis of NRT10, using both lensing 
and dynamics as independent probes of total mass, and 
providing a higher-redshift (z ~ 0.2) comparison to the 
z ~ galaxies previously studied 4 . 

2. DATA SAMPLE AND ANALYSIS 

2.1. Galaxy Samples 

The lensing galaxy sample is taken from the SLACS 
survey (A+09) 5 , which has been extensively analyzed in 
other works (e.g., Bolton et al. 2006, 2008; Gavazzi 
et al. 2007; Cardone et al. 2009; Treu et al. 2010; 
Grillo et al. 2009; Cardone & Tortora 2010). The lens 
galaxy redshift (zi) range is 0.05 < z\ < 0.5, with a 
median of zi ~ 0.2. Our sample is selected: 1) to have 
a measured Einstein radius Re 6 in Table 3 of A+09, 2) 

4 In the paper, we use a cosmological model with (O m , S1a> h) = 
(0.3, 0.7, 0.7), where h = ifo/lOOkms" 1 Mpc" 1 (Spergel et al. 
2007), corresponding to a universe age today of t un i v = 13.5 Gyr. 

5 See also http://www.slacs.org/ 

6 The Einstein radii are derived by fitting a singular isothermal 
ellipsoid (SIE) profile and are quoted adopting an intermediate- 
axis normalization. Five of the galaxies without a measured R^ 
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to be classified as elliptical or SO, 3) to have a measured 
V-band effective radius (measured at the intermediate 
axis). Of the 85 lenses from this latest SLACS release, 
66 passed our selection criteria. 

As a z ~ comparison sample, we use the collection of 
330 ETGs over the same mass range analyzed in T+09 
and NRT10. 

2.2. Stellar population analysis 

To estimate stellar mass-to-light ratios (T*) and star 
formation histories, we analyze spectral energy distri- 
butions (SEDs) based on broad-band Sloan Digital Sky 
Survey (SDSS) photometry (namely, ugriz). Our gen- 
eral procedure is to adopt a set of synthetic spectra from 
the prescription of Bruzual & Chariot (2003), a uniform 
metallicity Z, an age t characterizing the time of star 
formation onset, and an exponentially declining star for- 
mation rate with timescale r. For each galaxy, Z , t, and 
r are fitted parameters, with the determined Y* based 
on a Kroupa (2001) IMF; uncertainties on the estimated 
parameters have been quantified via Monte Carlo simu- 
lations: further details are provided in T+09 and NRT10 
along with explorations of systematic uncertainties and 
degeneracies. The only change here is to shift the spec- 
tral responses of the SDSS filters to correspond to the 
lens redshifts before convolving with the model SEDs. 
We have checked that imposing restrictions on Z or r, 
or adopting the stellar populations results from A+09 
or Grillo et al. (2009), does not qualitatively affect the 
results described below. 

2.3. Total mass and dark matter content 

We derive the deprojected total mass from dynamics 
and lensing observables, and separate the DM from the 
stellar components using the stellar mass estimates dis- 
cussed in the previous section. We adopt the stellar 
effective radius i? g' as the fiducial reference point for 
mass comparisons. For the lens galaxies, i? e ff is mea- 
sured in the F-band, corresponding approximately to 
the rest-frame B-band used for the local galaxies. In 
both galaxy samples, the mass constraints are gener- 
ally based on measurements at smaller radii (~ 0.1i? c ff 
and ~ 0.5-Reff, respectively), and therefore some extrap- 
olation is required. For the total mass distribution we 
adopt a singular isothermal sphere (SIS) with density 
p(r) = Og IS /(27rGr 2 ) (e.g. Koopmans et al. 2006, 
Gavazzi et al. 2007, Koopmans et al. 2009), where ctsis 
is an unknown normalization to be determined by fitting 
the observables. For the stars, we adopt a constant-Y» 
mass profile based on the Hernquist (1990) model. 

To estimate dynamical masses we have used the SDSS 
stellar velocity dispersions ctsdss ) measured within a cir- 
cular aperture of i? ap — 1.5". Briefly, we have adopted 
the spherical Jeans equation to derive the surface bright- 
ness weighted velocity dispersion <7 aP] sis within i? ap (see 
T+09 for further details), to be matched to ctsdss- As 
discussed in T+09, there is some degree of systematic 
uncertainty from assumptions of sphericity and orbital 
isotropy, which we can now check in the case of the lenses 
by using the independent lensing-based masses (which do 

have a nearby companion while, for other systems, the HST data 
do not have sufficient sensitivity to adequately perform the lensing 
model. 



have their own uncertainties from mass-sheet degenera- 
cies). 

For the lensing mass estimates, we have used the Ein- 
stein radius Re, to derive a model independent measure- 
ment of projected mass (M pro j) within Re, since Me = 

M proj (i? E ) = 7ri? E 2 S cri t, where E crit = c 2 D s /4ttGD 1 D 1s , 
with D s , Di and Di s the observer - source, observer - lens 
and lens - source comoving angular diameter distances, 
respectively. Finally we match the prediction of the SIS 
model projected mass, M pro jsis with Me to have a fur- 
ther constraint on the only free model parameter for each 
galaxy, ctsis- 

The best fitted ctsis can be derived independently us- 
ing either technique to estimate the best 3D deprojected 
mass profile which we extrapolate to r — R e s to ob- 
tain our reference mass values. Using this approach, we 
find that lensing and dynamics provide consistent results, 
modulo a ~ 10% (and a scatter of ~ 25%) higher mass 
from dynamics, corresponding to a change of 0.03 ± 0.10 
in /dm- Thus, we adopt a combination of the constraints 
for our final masses, by minimizing with respect to ctsis 
a combined x 2 function including one term for dynamics 
and one for lensing observables, given by 



2 _ / Q"qp — ctsdss \ ' 
X = ( S d ) 



Me — Mp ro j^sis 



(1) 



where 5d and Si are the uncertainties on ctsdss and Me, 
respectively 7 . 

In the following we will focus on the central 3D depro- 
jected /dm and the mean DM density within i? e ff) defined 
as (/9dm) = M D M/(4/37ri? cff 3 ) where M DM = M tot - M* 
at R e g is the DM mass. 

2.4. Cosmological models 

As in T+09 and NRT10, to interpret the observational 
results, we construct a scries of toy mass models based 
on ACDM cosmological simulations. For each bin in M*, 
we use the average i? c fj-age relations from the combined 
lens+local sample, and parameterize the virial DM mass 
by a star formation efficiency 6sf = M*/(f2bar-Mtot), 
where f^bar = 0.17 (Spergel et al. 2007) is the baryon 
density parameter. The halo densities are initially char- 
acterized as Navarro et al. (1997) profiles following an av- 
erage mass-concentration relation, adjusted by (1 + z)^ 1 
for the lens galaxies. A recipe for AC from baryon set- 
tling is then applied (Gnedin et al. 2004). The toy 
models for £sf = 0.03, 0.1, 0.3 are shown in both Figs. 1 
and 3. 

3. RESULTS: CORRELATIONS WITH SIZE AND 
FORMATION EPOCH 

In order to marginalize any correlations with M*, wc 
group the galaxies from both samples into bins of com- 
mon median mass: log M^/Mq ~ 11.6, 11.3 and 10.9. 
For the lens galaxies, the corresponding median redshifts 
are z me d = 0.28, 0.18, 0.13. The local galaxies sample ex- 
tends to even lower masses (log M*/Mq ~ 10.4) with no 
lens counterparts. 

7 Note that 8^ is given in Table 3 of A+09 and ranges from 2% 
to 19% (with mean 6%), while we have assumed a nominal 5% 
uncertainty on _Re which corresponds to a relative error of 10% on 
Me- However, the results are qualitatively unchanged if we would 
assume Sj = <5; . 
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Fig. 1. — Dark matter fraction within an effective radius (i? c ff) 
as a function of R s g . The lens and local galaxies are shown as filled 
and open symbols, respectively. For the latter, open symbols with 
error bars show the median and ±25% values. Typical ler uncer- 
tainties for individual galaxies are shown to the left. The differently 
colored symbols denote different bins of stellar mass, as labeled in 
the legend in the top. The second panel from the top shows the 
combined bins for the lens galaxies only and includes toy-model 
ACDM predictions: solid, long-dashed, and short-dashed curves 
show star formation efficiencies of egF = 0.03, 0.1, 0.3, respectively. 

Our first result (which we do not show for the sake of 
space) is that for the lenses, /dm increases on average 
with M* (see Cardone et al. 2009 and Cardone & Tor- 
tora 2010 for details). This result confirms that as in 
local galaxies (T+09), /dm is also a main driver of the 
fundamental plane tilt at z ~ 0.2. 

Next, following NRT10, we focus on correlations of the 
DM metrics with galaxy size and age. For the latter we 
adopt the look-back time to the formation epoch in order 
to put all the galaxies with different observed redshifts 
on a common reference frame. 

Fig. 1 demonstrates that there is a strong positive 
correlation between /dm and R e g, once the galaxies are 
divided into mass bins. This may be understood as a 
larger R e g enclosing a bigger portion of the DM halo; 
this "aperture effect" appears to be more dominant than 
the /dm correlation with M* . The local and lens samples 
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Fig. 2. — Mean DM density within 1 R e g as a function of R c g. 
See Fig. 1 for the meaning of symbols. The grey region shows an 
average of spiral galaxy results from McGaugh et al. (2007). 

appear reasonably similar, although the lens galaxies in 
the lowest mass bin are systematically higher, which is 
an issue we will discuss below. Both samples are in rough 
agreement with our ACDM toy model predictions (top 
panel) . 

Fig. 2 shows that (/9dm) strongly anti-correlates with 
-Roff- Again considering the aperture effect and assum- 
ing DM halo homogeneity, the implication is that we are 
measuring a mean DM density profile with radius, with 
a best fitted log slope of ~ —1.7. As discussed in NRT10, 
this steep slope is indicative of cuspy halos, perhaps as in- 
duced by AC. One could suspect that we are getting out 
what we are putting in, since our default galaxy model 
assumes an isothermal total density profile (with slope 
~ —2) in order to extrapolate measurements to r = R e g, 
but we have shown in NRT10 that the use of an alterna- 
tive constant-M/L profile yields similar results (modulo 
a difference of 0.1 — 0.2 in the slope), still fully consistent 
with a cuspy contracted halo 8 . 

In Fig. 2, we also see that the ETGs have DM densi- 
ties substantially larger than those of local spiral galax- 
ies, which have been suggested to follow a unified halo 
trend with dwarf spheroidals (Donato et al. 2009; Walker 
et al. 2010). This dichotomy is qualitatively consistent 
with other findings (Gerhard et al. 2001; Thomas et al. 
2009; NRT10; Cardone & Tortora 2010), and may imply 
different formation mechanisms. The difference might be 
assumed as simply caused by ETGs forming from denser 
late-type galaxies at earlier epochs, which would yield 
the corollary prediction that ETGs with younger stellar 
ages have less dense halos because their late-type pro- 
genitors were less dense. However, as we will see below, 
the opposite trend appears to be observed. 

Finally we consider the /DM-age dependencies in Fig. 3, 
again using separate M+ bins. The lens and local galax- 
ies match up remarkably well in general, showing a clear 
trend for lower /dm at older ages. The low-mass lens 

8 Note that the two massive galaxies J0157— 0056 and 
J0330— 0020 with logi? e gf ~ 0.9, which depart from the mean 
trend, are fitted by a possibly unrealistic supersolar metallicity; 
setting Z to the solar value, the galaxies' densities are incremented 
by 0.15 and 0.25 dex, respectively. 
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Fig. 3. — DM fraction within 1 R c ff as a function of galaxy stellar 
"age" . See Figs. 1 and 2 for the meaning of symbols. 

galaxies are predominantly young, which is probably a 
selection effect on apparent magnitude. The /DM _ age 
anti-correlation then causes the overall /dm for this mass 
bin to be high. In summary, at a fixed galaxy age and 
mass, the higher- z sample does not show any significant 
difference with the local galaxies. Since we have not ap- 
plied any evolutionary corrections, the implication is that 
in the last ~ 2.5 Gyr (on average) the galaxy populations 
have experienced no measurable evolution except stellar 
aging 9 . 

We also include the ACDM toy models in Fig. 3, and 
see that some /DM-age anticorrelation is expected, which 
can be traced to the anti-correlation between i? e g and 
age. However, there are indications in every mass bin 
that the observed /bM-age slope is steeper than in the 
models. Some systematic link between age and £sf is 
possible, but the models shown in Fig. 3 suggest that this 
would not be a strong enough effect, as changes in virial 
mass do not propagate strongly to changes in central 
DM content (and in fact earlier collapsing halos should 
have denser centers, which goes in the wrong way to ex- 
plain the observations). The alternatives as discussed in 
NRT10 are that AC is more effective in younger galaxies, 
or that older galaxies have "lighter" IMFs (e.g., Kroupa 
versus Salpeter 1955 for the younger galaxies). 

4. CONCLUSIONS 

We have analyzed the central DM content of a sam- 
ple of intermediate- z lenses from the latest release of the 
SLACS survey (A+09). Following the phenomenologi- 
cal framework introduced in T+09 and NRT10 we have 

9 Here the stellar mass bins can be biased at different redshifts 
due to the stellar population evolution and introduce a spurious tilt 
in the /DM _a g e relation. We have checked from the stellar models 
that the change in stellar mass due to stellar evolution might be 
not larger than ~ —0.01 in logM* (due to mass loss), which makes 
our stellar bins at different z fairly homogeneous. 



discussed scaling relations between DM fraction, galaxy 
size, and formation epoch. Gravitational lensing and dy- 
namical analyses are used to constrain the total mass 
profile, while synthetic spectral populations are used to 
infer the stellar mass and other stellar properties such as 
galaxy age. 

ETGs at z ~ 0.2 are found to be similar to local ones; 
future work will include extending the baseline to higher 
redshifts. The somewhat surprising findings of NRT10 
are now confirmed with an independent and arguably 
more robust data set. The DM fraction within R e g 
is found to strongly correlate with R e g, because larger 
length-scales probe a more DM dominated region. On 
these scales, the DM mean density decreases with R^g 
as (/0dm) oc -R c ff _ ' , which argues for cuspy DM halos 
for ETGs out toz~ 0.5 (as was discussed in NRT10 for 
local galaxies). At a fixed stellar mass and length-scale, 
we have found that the DM halos of ETGs are denser 
than those of local spiral galaxies, providing a critical 
test for the merging formation scenario (see also Car- 
done & Tortora 2010). 

Finally, we have confirmed our earlier finding that cen- 
tral DM content anti-correlates with stellar age. The 
strength of this correlation appears to exceed what is ex- 
pected from size-age effects. A fundamental connection 
between galactic structure and star formation history is 
implied, which we propose is a consequence of variations 
with formation epoch of either DM halo contraction or 
stellar IMF. 

In future work, we plan to investigate the impact on 
these results of more complex total and DM galaxy pro- 
files (e.g. Cardone et al. 2005, Tortora et al. 2007) along 
the lines of recent work in Cardone et al. (2009) and 
Cardone & Tortora (2010). New high-quality data are 
also expected with the advent of future surveys both in 
the local Universe and at larger redshifts. Such surveys 
will include larger samples of gravitational lenses along 
with more detailed spectroscopic information, and could 
be used to verify and extend the results presented here, 
providing a clearer picture of the physical processes of 
ETG assembly. 

We thank the anonymous referee for the suggestions 
which helped to improve the paper. CT was supported 
by the Swiss National Science Foundation. AJR was 
supported by National Science Foundation Grants AST- 
0808099 and AST-0909237. 

Note added in proof. During the final phase of this 
manuscript, there appeared closely related work from 
Auger et al. (2010a, b, A+10a and A+10b respectively, 
hereafter), using a combination of strong lensing, and 
stellar dynamics and populations, to analyze virtually 
the same lens-galaxy sample. A key point of agreement 
is that A+10a found the strongest non-trivial correlate 
with central /dm is R c e, which we think demonstrates 
that size rather than mass variation is the main driver 
for the fundamental plane tilt. 

A+lOb also constructed ACDM models, adding con- 
straints from weak lensing to break the degeneracies be- 
tween £sf, IMF, and halo model (AC or no-AC). Given 
the same IMF and halo assumptions, they found SF 
values that are on average higher than both ours (e.g. 
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6sf ~ 0.3 — 0.4 versus ~ 0.1 for Salpeter+no-AC), and 
typical values found in other studies. Although these 
differences may seem large, they involve relatively small 
changes in the central DM properties which could be 
driven by differences in the modeling methods. We con- 
sidered three-dimensional DM masses within i? e ff, while 
A+10b analyzed projected masses within R c s_/2 in com- 
bination with large-radius constraints from weak lensing, 
so a direct comparison is not straightforward. We hope 
to track down the reasons for these discrepant conclu- 
sions in future work. 

A+lOb found that a Salpeter+no-AC model is pre- 
ferred over a Chabrier IMF (with or without AC) using 
direct model fits. Without weak lensing constraints, our 



modeling permits either Salpeter+no-AC or Kroupa+AC 
solutions (see NRT10 for further details), while our Fig- 
ure 2 (or Figure 9 in NRT10) does provide indirect evi- 
dence for very cuspy halos with AC (and a Kroupa IMF). 
They also found that the IMF may become heavier with 
galaxy mass (77 > 1 in their notation) . We find that stel- 
lar age is the more important parameter in this context, 
but if an IMF-mass relation is demanded, then 77 < 1. 
This difference appears to be caused at least partially by 
different stellar populations models: although we found 
that varying these models did not alter the basic age 
trends, the mass trends are weaker and more sensitive to 
the models. 
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